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Abstract 



We show, essentially by theoretical means, that for a site with the 
chosen simple geometry and mechanical properties (horizontal, homo- 
geneous, soft viscoelastic layer of infinite lateral extent overlying, and 
in welded contact with, a homogeneous, hard elastic substratum of half- 
infinite radial extent, shear-horizontal motion): 1) coupling to Love 
modes is all the weaker the farther the seismic source (modeled as a 
line, assumed to lie in the substratum) is from the lower boundary 
of the soft layer, 2) for a line source close to the lower boundary of 
the soft layer, the ground response is characterized by possible beating 
phenomena, and is of significantly-longer duration than for excitation 
by cylindrical waves radiated by deep sources. Numerical applications 
of the theory show, for instance, that a line source, located 40m below 
the lower boundary of a 60m thick soft layer in a hypothetical Mexico 
City-like site, radiating a SH pulse of 4s duration, produces substan- 
tial ground motion during 200s, with marked beating, at an epicentral 
distance of 3km. This response is in some respects similar to that 
observed in real cities located at soft-soil sites so that the model em- 
ployed herein may help to establish the causes and pinpoint the major 
contributing factors of the devastating effects of earthquakes in such 
cities. 

Keywords: site response, regional path effects, source position. Love modes, 
interference maxima, duration, beatings. 
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1 Introduction 



This investigation is relevant to several topics of broad interest in seismic 
wave propagation: 

(a) regional path effects in connection with seismic response in urban 
environments [l], [S], |H], 0, 0, 0, 0, 0, 0, M 

(b) effects of the underlying soil heterogeneities, lateral variations of the 
underlying soil layer, and built environment on seismic ground response at 
various (particularly urban) sites [11], [E], [E], 0, [H], [E], [E], [T7]-F^- 

0, m, m 

(c) analysis of surface wave response on the ground to determine the 
structure and composition of the crust and underground fault zones 

m, m 

(d) analysis of surface wave response on the ground to identify earthquake 
sources [22], |23| . 

Research on topic (a) was rekindled by efforts to explain some puzzling 
features of the devastating Michoacan earthquake which struck Mexico City 
in 1985. Other than the fact that the response in downtown Mexico varied 
considerably in a spatial sense, was quite intense and of very-long duration 
|24| . [2H] (as much as ~3 min) at certain locations, and often took the form 
of a quasi-monochromatic signal with beating, a remarkable feature of this 
earthquake was that such strong (in the sense just mentioned) response could 
be caused by a seismic disturbance so far (its epicenter was in the subduc- 
tion zone off the Pacific coast approximately 350 km) from the city [2^], |27j . 
0, Pn]- A part of the cause of the large intensity and long-duration was 
attributed in to multipathing between the source and the site. This hy- 
pothesis was further explored in [2S], [2Zj while being associated with surface 
wave propagation of the Rayleigh and Love types, presumably between the 
source and the entry to the Mexico City basin, via the intervening crust. 

In a rather complete (other than the neglect of attenuation) 3D numer- 
ical study 0, the long duration and large amplitude of response at various 
distances from subduction zone earthquakes in Mexico were attributed to 
the entrapment of the seismic disturbance in an acccretionary prism (wedge- 
shaped heterogeneity) of the crust and its subsequent propagation to the 
point of observation. The authors of this work later 0, stressed the role 
of higher-order surface waves which propagate in the relatively-high Q layer 
of the Trans Mexican Volcanic Belt (TMVB) underlying the soft clay basin of 
Mexico City in addition to that of the accretionary prism in producing large 
response (particularly with respect to duration) in the city. More recently, 
an analysis ^ of seismograms recorded at various sites in central Mexico, for 
earthquake sources located in the subduction zone off the Pacific coast, have 
shown that the crustal structure (including that of the TMVB, composed 
of low-velocity volcanic lava and tuff overlying higher velocity limestone) 
between the source and observation points "acts as a waveguide for surface 
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waves coming from distances greater than 200km", leading, by an unex- 
plained mechanism, to amplification and increase of duration of motion at 
various sites, this being thought to account for at least part of the anomalous 
response in Mexico City to remote seismic disturbances. Numerical results 
obtained in earlier studies (e.g., (Hj, ^01) with a rather complete, 2D hy- 
brid model of the propagation path between the source and the Mexico-City 
basin, and of the action of the basin on the incident wave, also stressed the 
important role of regional path effects on anomalous response. 

Anomalous response in other cities such as Beijing, Bucharest, Rome, etc. 
has been studied in great detail, principally in numerical manner, within the 
framework of the UNESCO-IGCP project 414 [12], [H]. The features of this 
response were attributed to the specifics of the source parameters, regional 
path effects, and the specifics of the soil distribution and geometry in the 
urban basins (see next paragraph). These findings have been substantiated 
in a more recent study [TH] . 

Topic (b) deals with a class of alternative or complementary (so-called 
local) paradigms for explaining seismic motion in urban sites built on soft 
soil. Even though the anomalous response in 1985 in Mexico City originated 
in a subduction zone source whose epicentral distance was some 350 km from 
the city, it has been common to seek explanations of this response (and others 
such as in Nice pH], ^H] and Bucharest |T3j) by employing models involving 
vertically-propagating or nearly vertically-propagating plane waves. This 
requires that the focal distance of the source to the surface be large and that 
the epicentral distance from the source to the city be rather small. Although 
both of these conditions are often not met in practice (and, in particular, as 
concerns the 1985 Michoacan earthquake), the vertically-propagating plane 
wave sollicitation usually prevails in the theoretical/numerical studies [Tlj . 
HI; [2Hj; apparently because it simplifies the analysis (another reason is that 
it facilitates comparison with the so-called ID model of normally-incident 
plane waves on a vertically-layered half space) . This has the effect of putting 
the focus on what occurs in the structure vertically below the city, namely, 
on the soft basin on which the earthquake-prone cities are built. Thus, a 
considerable amount of studies (see jSHI; 1211; HE] for comprehensive reviews) 
examine the (local) effect of the soft basin on the incident wave, but at 
present, it is thought that local effects account for only part of the anomalous 
response (2^], (22j, plj, |H], ^H|' Another idea that has been explored in 
the past few years is that the buildings of the city, in interaction with the 
soft soil and with each other, may also amplify and lengthen the duration 
of the ground motion (see ^21; ^Hj for reviews of this subject). All of these 
studies (including or excluding the buildings) point to the central role of 
surface waves, qualified either as locally-launched surface (e.g.. Love) waves 
(at the basin edges or at heterogeneities of the soft soil) [20], [20], [0], HE] 
or as quasi-Love waves (excited at the base of the buildings and re-amplified 
by interaction with neighboring buildings ^2], [H^, QH]) as a possible causal 
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agent of anomalous response, but little jHS], if any, theoretical evidence has 
been given to back up these assertions. 

Topic (c) is classical in seismological geophysics [SS]. The seismic signals 
associated with various types of surface (e.g., Love and Rayleigh) waves 
are oft-used tools for reconstructing features of the earth's crust such as 
thickness, composition (e.g., vertical layering characteristics [33], [31]) and 
even lateral heterogeneities [32], More recently [20], [2^, it has been 

shown that seismic sources in the neighborhood of fault zones (FZ, i.e., soft 
nearly-vertical layers surrounded by relatively hard soil) excite surface waves 
(qualified as "trapped") in the vicinity of the FZ which propagate to the 
ground where they can be detected and used to furnish information on the 
physical and geometrical characteristics of the FZ. To treat these inverse 
problems in a fully unambiguous manner requires a thorough understanding 
of the way the seismic source interacts (notably how accurately one must 
know the position and characteristics of the source) with the inhomogeneities. 

Topic (d) is also a classical one in seismology, the main concern being 
to localize and qualify (e.g., determine the moment tensor of) earthquake 
sources [22], [23]- As the seismic wave, including its surface- wave compo- 
nents, travels laterally (sometimes over long distances) in and along the 
crust before reaching the measurement locations on the ground, the inverse 
problem is difficult to solve if the crustal features (which can include lat- 
eral heterogeneity) are not known beforehand. In any case, it is important 
to determine the influence of errors of the crustal model on the reconstruc- 
tion of the source location and moment tensor, and to do this requires an 
appropriate theoretical analysis. 

The theoretical-numerical investigation herein is focused on topics (a) 
and (b). In contrast to the inverse-scattering topics (c) and (d) (to which 
our analysis could be applied) wherein the response is known and the prop- 
agation medium and/or the source are to be determined, the problem we 
are faced with herein deals with forward- scattering: given the seismic source 
and the characteristics of the propagation medium, determine the response 
(displacement in the frequency and/or time domain) on the ground. More 
specifically, we shall be concerned with a (deceivingly-) simple canonical 
scattering problem: that of a cylindrical SH pulse wave impinging on a soft 
homogeneous layer, the latter being horizontal, of infinite lateral extent, 
bounded above by the ground and below by an interface with a half space 
filled with hard homogeneous rock. The questions we address, and that we 
think can be answered with the help of such a simple model, are: 

(i) is it possible to obtain anomalous (in the sense mentioned above in 
connection with the Michoacan earthquake) response without any lateral 
heterogeneity (arising from volumic inclusions or uneveness of interfaces) in 
the underground medium? 

(ii) what is the relation of ID to 2D response and how adequate is it to 
model the general response of the configuration by its response to a (nearly) 
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vertically-incident plane wave? 

(iii) how does the focal distance of the source affect the response? 

(iv) how does the epicentral distance affect the response? 

(v) how does the contrast of mechanical properties between the layer and 
the half space affect the ground response? 

(vi) how does the thickness of the layer affect the response? 

(vii) how do the spectral characteristics of the incident pulse affect the 
response? 

It will be shown that a source radiating cylindrical waves in a fully-elastic 
soft layer/hard half space medium produces a ground response which is the 
sum of three terms corresponding to various combinations of two types of 
waves in the soft layer (SL) and hard half space (HHS): 

(1) standing body waves (SBW) in the SL and body waves (BW) in the 
HHS, 

(2) standing body waves in the SL and surface waves (SW) in the HHS, 

(3) standing surface waves (SSW) in the SL and surface waves in the 
HHS. 

Only type (2) waves correspond to Love modes (at the resonance frequen- 
cies of these modes) and the conditions for optimal excitation and maximal 
contribution of these modes will be rendered explicit. It will be shown that 
large-duration (i.e., anomalous) response generally requires a preponderant 
contribution of at least one (usually the lowest-order) of the Love modes to 
the overall response. The type (1) waves dominate in the situation in which 
the focal distance is large and do not usually produce long-duration response, 
although they can produce strong (although normal) response when the con- 
trast of mechanical properties between the SL and HHS is large. Beating 
phenomena will be shown to be a consequence of interference between type 
(1) and type (2) waves which both lead to maxima in response at nearly 
the same (low) frequency. Type (3) waves turn out to have negligible con- 
tribution to overall response. Most of these features carry over to the case 
in which the layer is lossy. The practical consequences of these results, in 
relation with topics (a) and (b), will be discussed. 

2 Description of the configuration 

Fig. 12.11 represents a cross-section (sagittal plane) view of the site. Tg 
is the ground, assumed to be flat and horizontal, above which is located the 
air medium, assumed to be the vacumn. is the laterally-infinite domain 
occupied by the mechanically-soft layer and h is its thickness. Qq is the semi- 
infinite domain (substratum) occupied by a mechanically-hard medium, and 
Th the flat, horizontal interface between the layer and the substratum. A 
0x1X2X3 cartesian coordinate system is attached to this configuration such 
that O is on the ground, X2 increases with depth and X3 is perpendicular to 
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the (sagittal) plane of the figure. With ij the unit vector along the positive 
Xj axis, we note that the unit vectors normal to Tg and are 12- The media 
filling Qq and are and respectively and the latter are assumed 
to be initially stress-free, linear, isotropic and homogeneous. We assume 
that M° is non-dissipative and is generally (unless specified otherwise) 
dissipative. 

The seismic disturbance is delivered to the site in the form of a shear- 
horizontal (SH) cylindrical pulse wave radiated by a line source (perpendic- 
ular to the sagittal plane) located at x" := (xf,x|), with, by hypothesis, 
X2 > h (i.e., x'^ G flo). The SH nature of this wave means that the motion 
associated with it is strictly transverse (i.e., in the 3:3 direction and indepen- 
dent of the X3 coordinate). Both the SH polarization and the invariance of 
the incident wave with respect to X3 are communicated to the fields that are 
generated at the site in response to the incident wave. Thus, our analysis 
will deal only with the propagation of 2D SH waves (i.e., waves that depend 
exclusively on the two cartesian coordinates xi, X2 and that are associated 
with motion in the X3 direction only). 

We shall be concerned with a description of the elastodynamic wave- 
field on the ground (i.e., on Tg) resulting from the cylindrical seismic wave 
sollicitation of the site. 



3 Governing equations 

3.1 Space-time domain wave equations 

In a generally-inhomogeneous, isotropic elastic or viscoelastic medium M 
occupying M'^, the time-domain wave equation for SH waves is: 

V • (^(x, u)Vu{jc, t)) - p(x)a2tx(x, t) = -p(x)/(x, t) , (3.1) 

wherein u is the displacement component in the i3 direction, / the component 
of applied force density in the 13 direction, /i the Lame descriptor of rigidity, 
p the mass density, t the time variable, u; the angular frequency, the n— th 
partial derivative with respect to t, and x = (xi, X2). Since our configuration 
involves two homogeneous media and the applied force is assumed to be non- 
vanishing only in Oq, we have 

(c"(w))'V-V^x™(x,t)-a2tx-(x,t) = -/(x,t)5^o ; xgQ^, (3.2) 

wherein m superscripts designate the medium (0 for or 1 for M^), 6mO = 
1 for m = and equal to zero otherwise, and c™ is the generally-complex 
velocity of shear body waves in M™, related to the density and rigidity by 

{c-{u.)f = ^ , (3.3) 

it being understood that p'^,pr^{uj) ; jn = 0, 1 are constants with respect to 
x. 
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3.2 Space-time domain representation of the impulsive force 

In all that follows we choose the pseudo-Ricker type of impulse line source 
function 



/(x, t) = -5(x - x^)3^ -2a2 (i - 2a^ {t - tof) e-"'(*-*°)' 



d_ 
'dt 

5(x - x^)12a^ [-3 {t - to) + 2a^ {t - tof 



-aHt-tof 



, (3.4) 



wherein a = ttuq and Iq = and 5{ ) the Dirac delta distribution. 

3.3 Space-frequency domain wave equations 

The frequency-domain versions of the wave equations are obtained by 
expanding the force density and displacement in Fourier integrals: 



/(x,i) 



(3.5) 



so as to give rise to the Helmholtz equations 



V • Vt/'"(x, uj) + ik"'ioj)f n™(x, a;) 



f{x,uj)6mO ; Vx G ; m = 0, 1, (3.6) 



wherein 



UJ 



(3.7) 



is the generally-complex wavenumber in M'". Actually, due to the assump- 
tions made in sect. E 



,0/ N IP 



(3. 



(i.e., is real). 



r(x,u;) = S(cu)5(x-x^) 



(3.9) 



wherein S{io) is the spectrum of the incident pulse. In fact, the spectrum 
corresponding to the chosen (see (j3.4|l pseudo-Ricker impulsive force is 



(3.10) 
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3.4 Material constants in a dissipative medium 



A word is now in order about the dissipative nature of the layer. When a 
medium M is lossy, the wavenumber therein is complex and can be written 
(omitting, for the moment, the lo dependence) as 



where, by convention. 



Rk = k' >0 , = /c" > . 



(3.11) 
(3.12) 



We now refer to (|3.7|l and note that complex k implies complex /i, due to the 
fact that it is advisable to consider the mass density to be a real quantity. 
Thus, we write 

/i = ^' -i/i" . (3.13) 

In order to retain the positive real aspect of the rigidity for elastic materials, 
we take 

^l_i = >0 , (3.14) 

and inquire as to the sign of the imaginary part of fi. Introducing (|3.13() into 
H3.7P gives 





-1/2 















-1/2 



(3.15) 



We assume, as is generally the case for moderately-dissipative media, that 
|^"/a*'| << 1) so that a Taylor series expansion of [ J"^/^ limited to the first 
two terms yields 



k = k' + ik'' 



wherein, by definition. 



P 



l + i 



1/2 



Making use of (|3.12p and H3.13|) thus necessarily leads to 

= < 0. . 



(3.16) 
(3.17) 
(3.18) 



We define the positive real quantity known as the quality factor Q by the 
ratio 



Q 



(3.19) 



and note that it is infinite for a lossless medium such as (because /i" = 
in this case). Furthermore, the complex wavenumber becomes 



(3.20) 
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from which we find 

A question arises as to the proper definition of the complex body wave ve- 
locity c in M. We write 

c = c' - ic" , (3.22) 

and require 

^c = c' >0 , (3.23) 

due to the fact that the body wave velocity is positive in a non-lossy medium. 
We have 

,/ ^ ^ ujc' + iujc" 

k = k' + ifc" = - = — = ^ , 3.24 

c c — ic |c|^ 

from which we see that in order for = k" > 0, we must have 

= -c" < . (3.25) 

The remaining question is that of the w-dependence of fj, and Q (the w- 
dependence of k and c follows from that of n and Q). In seismological 
applications involving viscoelastic media the quality factor is found to be ei- 
ther constant or a weakly- varying function of frequency We shall assume 
that Q^{lli) = =const., and it can be shown [SEl that this implies 



/iH^) = /^re/ ( ^ ) (3-26) 

wherein: oj^ef is a reference angular frequency, chosen herein to be equal to 
9 X lO^^Hz. Hence 



"ref 

and 



cHu.)=clf{^]- -\ (3.27) 



Note should be taken of the fact that even though is non-dispersive (i.e., 
does not depend on uj) under the present assumption, the phase velocity 
is dispersive. 

3.5 Boundary and radiation conditions 

We assume the two media to be in welded contact so that the displacement 
and the normal components of stress are continuous across the interface F/j: 

ui(x,a;) - u°(x,cj) ; X E T/j , (3.29) 
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/i^(u;)9„uXx,u;) - (w)a„n"(x, cj) ; x G r,^. (3.30) 

Since the air/layer interface Tg (i.e., the ground) is assumed to separate the 
vacumn from an elastic medium, the normal component of stress must vanish 
on this boundary, i.e., 

fi^{uj)dnu'^{^,uj) = ; X G , (3.31) 

wherein 9^ = i2 • V = ■ The uniqueness of the solution to the forward- 
scattering problem is assured by the radiation condition in the substratum: 

u^{x,uj) ~ outgoing waves , ||x|| — > oo, X2> h . (3.32) 

3.6 Statement of the boundary- value (forward-scattering) prob- 
lem 

The problem is to determine the time record of the ground displacement 
field u^{xg,t) (with x^ := (xi, 0)) from the spectrum of the ground displace- 
ment u^{^g,uj) via the Fourier transform 

/oo 
u^{yLg,uj)e-''^^dio . (3.33) 
-oo 

Note that due to the fact that u^{xg,t) is a real function, we must have 

[u^(xg,u;)]* = u^(xg, -cj) , (3.34) 

(wherein the symbol * designates the complex conjugate operator) from 
which it follows that 

POO 

u^(xg,t) = 23f? / u^{xg,uj)e-''^^dio . (3.35) 
Jo 

4 Exact solutions in the frequency domain by sep- 
aration of variables 

4.1 Preliminaries 

Although the material in this section (4) is classical as regards the way 
of obtaining plane wave integral representations of the fields, the way these 
integrals are decomposed, analyzed, and computed is different than in previ- 
ous investigations (e.g., [HZ|,|HHI, 1121,1101,1111, |l2l,|lHI,|llI,|lSI,|lSl, |1ZI,|1HI, 
0ni , ISni , |S2 ; the considerable quantity and variety of these publications at- 
tests to the richness and importance of the subject, and to the fact that 
certain features of the latter certainly remain to be discovered). In a first 
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subclass of these investigations, the plane wave integrals (with the horizon- 
tal wavenumber as the variable of integration) are reduced to residue series 
(so-called modal series) plus branch cut integrals which are usually neglected 
if the source-to-observation point is large compared to the wavelength. In a 
second subclass of the aforementioned investigations, various devices are em- 
ployed to evaluate in numerically-efficient, accurate, or asymptotic manner 
the plane wave integrals. Our contribution is essentially of the second vari- 
ety, but numerical efficiency (more important in the inverse problem context) 
is of less interest to us than the physical significance of the terms entering 
into our choice of the decomposition of the integrals. 

4.2 Frequency- domain solutions in the absence of the layer 
and the free surface 

In the absence of the layer and the free surface, the problem is that of 
the radiation of a SH wave from a line source in 2D free space (R^) occupied 
by the homogeneous medium M^. We term this radiated wave the 'incident 
wave' and designate it by n*. 

By applying separation of variables in the cartesian coordinate system to 
the Helmholtz equation and using the radiation condition, it can be shown 
that takes the form jH21 

471" J -00 '^2 

or 

nHx,a;) = ^5(^)<) (fcO(^)||x - x^||) , (4.2) 
wherein Hq^\ ) is the zeroth-order Hankel function of the first kind and: 

kl{uj) := ^{y{oj)f -kl , ^k{{oj) > Q , ^k{{uj) > Q , j = 0,l. (4.3) 

We shall make use in sect. 14.41 of the form taken by in the region : = 

{xl> X2> h \ Mxi e M}: 

/oo 
A0(A;i,u;)en'=i^i-'=2('^H^^^ . Vx G , (4.4) 
-CXD 

wherein 

A\ki,Lo) = S{lo) -i[fcix;-fcgMxi] _ (4 5) 

47rK2 [uJ) 

4.3 Field representations in cartesian coordinates for the 
configuration including the layer and the free surface 

When the layer and free surface are present, the incident field described 
in the previous section cannot proceed in unobstructed manner, i.e., it gives 
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rise to a 'diffracted' field (indicated by the superscript 'd') so that by re-use 
of separation of variables in cartesian coordinates and the radiation condition 
we are led to represent the total fields in the substrate and the layer by 

n°(x, u) = u*(x, u) + u°'^(x, cj) , (4.6) 

u^(x,cj) = ui^(x,cj) . (4.7) 

wherein: 

/oo 

5°(A;l,a;)eH'=l=^l+*^2H(^2-/^)J^^^ . > /i , Vxi G M, (4.8) 
-00 

n^'^(x,a;) = 

^i(fci,w)e^['=i^i-'=2('^)^2] +Bn'fci,'^)e'[*^''^''^''^^'^^'''^]) dki ; 

< X2 < /i , Vxi G M , (4.9) 

it being understood that the diffraction coefficients B^, , are, as yet, 
undetermined. 

4.4 Determination of the diffraction coefficients and frequency 
domain fields by application of the boundary conditions 

The free-surface boundary condition entails: 

^^(/ci,u;) = 5^(A:i,u;) , V/ci G M , (4.10) 

whence 

/oo 
A^{ki,uj)e'^^'''' cos{k\{oj)x2)dki ; < X2 < /i , Vxi G M. 
-00 

(4.11) 

The continuity of displacement condition leads to: 

B^{ki,uj) - 2A^{ki,uj)cos{k\{uj)h) = -A°(A;i, w)e-^^2('^)^ ; Vfci G M , 

(4.12) 

whereas the continuity of normal stress boundary condition implies: 

i/(w)A:^(a;)5°(A;i,a;) -h2^ifc^(u;)^^(/i:i,a;)sin(A;^(u;)/i) = 

i^x'^{uJ)kl{uJ)A\kl,uJ)e-'^'^^^'^^ ; Vfci G E . (4.13) 
The solution of this linear system of equations is: 

5°(/fci,w) = A0(fci,u;)e-^'=2H'^x 

' -jJ^ (^)fc2 (^) sin(fc2 + i/i°(c^)/c2 (^) cos(fc| {u)h) \ . ^ 

/i^(i^)A;2(w) sin(A;|(u;)/i) -|- i/i'^(Lij)/c2(<-o') cos(A;2(w)/i) / ' ' 

(4.14) 
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so that the solutions for the fields in the frequency domain are: 



; VA;i G M . (4.15) 



u''(x, U>) = S{U!) 



[fci(xi-zf)+fcOH|x2-x||]^^_^ 



+ 



Airk^iuj) 



if/'{uj)k2{uj) cos{k2{uj)h) — fi^{uj)k2{io) sm{kl{u;)h) 
i^^{ijj)k2(uj) cos{k\{ijj)h) + ^^{ijj)k\(uj) sm{k\{uj)h) 

gi[fci(xi-a;f)+fcOH(x2+:.|-2h)]^^^ ; Vx G , (4.16) 



u-^(x, w) = ti^'^(x, w) 



47r/c0(u;) 
cos 



2i^j^{uj)kl{uj) 



ifi^{uj)k2{uj) cos{kl{uj)h) + fi^{uj)kl{uj) sm{kl{uj)h) 



{kl{Lo)x2) e*['=i(^i-^i)-'=§H('^-^i)](iA;i ; Vx G f^i . (4.17) 
Finally, the frequency- domain ground response takes the form: 

S{UJ) 



2in^{uj)k^{io) 



A7rk^{uj) [ifi^{Lu)k^{Lu)cos{kl{Lu)h) + fi^{Lu)kl{Lu) s'm{kl{uj)h) 

gi[fcl(xi-x?)-fcOH(fc-xi)]^^^ . ^ _ 



5 Structure of the frequency-domain response in 
the case of a non-lossy layer 

5.1 Frequency domain response in the layer 

When the layer is free of dissipation, i.e., elastic, then /x^ is real and does 
not depend on w, and k^i^uj) is real (recall that we assumed the substratum 
to be elastic, which means that fi^ is real and does not depend on u>, and 
k^{u!) is real also). Consequently, in the integrals of the previous section 
we encounter intervals of ki over which and k^ are either purely real or 
purely imaginary: 



ki{uj) 



{y{u:)Y-ki 



■ \ki\ < k^{uj) ; (J > 



(5.1) 



k2{uj) = iK2{uj) := i 



kj - {kj{ioy) 



\ki\ > y{uj) ; w > , (5.2) 
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It is important to note that the terms 'soft layer' and (relatively) 'hard 
substratum' have the following meaning in the present context: 



so that (|4.17p can be expressed as: 

(x, oj) = ll (x, u>) + l\ (x, uj) + l\ (x, uj) ; Vx G ili 

with: 



(5.3) 
(5.4) 

(5.5) 



dti}(x,Xg, A:i,u;) = 

_ Flikuu:)cos {K^{co)x2) e*['=i(-i-i)-^§H('^-i)]dA:i , (5.6) 



Flikuoo) -- 



/|(x,u;) 



dul{x,Xg,ki,uj) = 



(5.7) 



+ / 



27r 



+ / 

fei Jko 



F2{ki,uj) cos [K2ii^)x2) X 

g[ifci(xi-xf)+«:0(a;)(h-x|)]^^^ ^ g-j 



F2^(/ci,u;) 



4(x,u;) 





/•oo 




^ Jk^ 


./ — oo 




S{u;) 


r 


27r 


J —oo 


\ — 





cos{Kl{uj)h) + ^i^Kl{uj)sm{Kl{uj)h) ' 



du\{'X-, Xg, fci, w) 



(5.9) 



+ 



F3^(xg, /ci, w) cosh (K2(a;)x2) x 



/i^Kgli^) cosh.{K\{uj)h) + /i^K2(w) sinh(K2(tij)/i) 



(5.11) 
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We write: 



dG\{^g,ki,u;) = -^F/(yfci,w)e-*['=^-1+^"(-)('^--'^)]dfei , (5.13) 

TT 

which, together with (|5.6p . express the fact that a part (i.e., /|) of the field 
in the layer is composed of a sum of standing body waves (SBW), each of 
which is the sum of two plane body waves having wavevectors with the same 
length. 

In the same manner, we write: 

dul{^,^g,ki,uj) = dG^(xg, A:i,w)e' [^■1^1+^2 (^H + 

(iG^(xg,A:i,w)e'[^'^^^~^2('")^2] , (5.14) 

dGl{x„,h,oj) = -^F^{ki,oj)e-'''''''^+^>^(^-'''^)dki , (5.15) 
vr 

which, together with (j5.8|l . express the fact that another part (i.e., I2 ) of the 
field in the layer is again composed of a sum of standing body waves, each 
of which is the sum of two plane body waves with wavevectors having the 
same length. Note however that neither the wavevectors nor the amplitudes 
of these SBW are the same as those of the SBW (henceforth termed SBWl) 
in /| (because the range of integration in the latter is different from that 
in I2). In fact, (|5.15jl tells us that the amplitudes dG^ of the SBW in I2 
(henceforth termed SBW2) decrease exponentially as the focal distance (i.e., 
X2) increases, so that the SBW2 make themselves felt all the less the farther 
the source is (in the vertical direction) from the ground. On the other hand, 
the amplitudes of the SBWl are sinusoidal functions of focal distance, so 
that the SBWl can possibly make themselves felt strongly for a large variety 
of source locations. 
Finally, we write: 

dul{^,^g,kuu;)= dGl (x, , ki , u;)e^^i^i +^^2 + 

(iGj(xg,A;i,w)e^'=^^i-^2('^)^2 , (5.16) 

dGl{xn,ki,oj) = -^Fi{ki,oj)e-''''''i+''>'>(''-''^^)dki , (5.17) 

TT 

which, together with (|5.inp . express the fact that the third part (i.e., /|) of 
the field in the layer is composed of a sum of standing surface waves (SSW), 
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each of which is the sum of two plane surface waves with wavevectors having 
the same length (note that each such plane surface wave is an inhomogeneous 
wave (with complex wavevector) whose phase is constant on xi =const. 
and whose amplitude either increases or decreases as X2 approaches some 
horizontal surface X2 =const.). Eq. H5.17|l tells us that the amplitudes dG^ 
of the SSW in Ig decrease exponentially as the focal distance increases, so 
that the SSW make themselves felt all the less the farther the source is (in 
the vertical direction) from the ground. 

The main conclusion of this discussion is that for focal distances of the 
source that are sufficiently large, the field in the layer is essentially given by 
Ii and is expressed by a sum of SBWl. This corresponds more or less to the 
situation in the quasi-lD analysis of the forward-scattering problem, but, as 
we shall see further on, it is, by no means, a valid picture of the response of 
the configuration when the focal distance of the source is not large. 

5.2 Frequency domain response in the hard half space 

We shall concentrate our attention exclusively on the diffracted field in 
the subdomain although the essence of what will be written applies to 
the whole half space • Proceeding as in sect 15. II we find: 



u"'^(x, io) = /^x, Lo) + /^(x, Lo) + /^(x, u) ; Vx G 



(5.18) 



with: 





K°{uj) ifi^K^iio) cos{Kl{uj)h) + p^Kl{uj) sin(Ki(w)/i) 



(5.20) 





F^{ki,uj)e^M^^-^\)-<^l^^){^^+''2--i'^)]dki , (5.21) 



1 -/i°K^(cj) cos{K^{uj)h) - iJ^Kl{uj) siu{Kl{uj)h) 
;^(cj) -pPi4{uj) cos{Kl{uj)h) + ii^Kl{uj) sin{Kl{uj)h) 



(5.22) 
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I3"(x,a;) 



47r 









J — oo 















+ 



dn3(x,Xg, A;i,u;) = 



1 — /i''K2(ti^) cosh(K2(u;)/i) + k\{uj) sm\i{K\{u))h) 
K2(w) — //'^«;2('^) cosh(K2(w)/i) — ii^k\{uj) sinh.[K\{uj)h) 



We write: 



(5.24) 
(5.25) 



dG?(x„,fci,a.) = ^F0(fci,u;)e^^['^^^i"^°('^K-i-2^)]dfci , (5.26) 
47r 

which, together with (|5.19|1 . express the fact that a part (i.e., Ii) of the 
diffracted field in the half space is composed of a sum of plane body waves 
(BW). Thus, to each horizontal wavenumber ki in the interval [—kP,k^], 
correspond a SBWl in and a BW in . 
In the same manner, we write: 



(5.27) 



dGO(x„A;i,u;) = ^F^{ki,Lo)e-''^'^'^-^°^<^'^-^'')dki , (5.28) 

which, together with (|5.2ip . express the fact that another part (i.e., I2 ) of the 
diffracted field in the half space is composed of a sum of plane surface waves 
(SW), henceforth denoted by SW2. Eq. (|5.28p tells us that the amplitudes 
dGg of the SW2 in decrease exponentially as the focal distance increases, 
so that the SW2 make themselves felt all the less the farther the source is (in 
the vertical direction) from the ground. On the other hand, the amplitudes 
of the BW in are sinusoidal functions of focal distance, so that these BW 
can make themselves felt strongly for a large variety of source locations. In 
addition, we note that to each horizontal wavenumber ki in the intervals 
[-k^, -k°] and k^], correspond a SBW2 in Qi and a SW2 in . 
Finally, we write: 

dul{^,^g,k,,oo) = dGl{^g,ki,oo)e'''^^^-'^»' , (5.29) 
dGlixg,ki,uj) = ^F3i(A:i,cj)e-*'=i^i-'^°('^K^i-2/^)dA:i , (5.30) 
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which, together with (|5.28p . express the fact that the third part (i.e., of 
the diffracted field in the substratum is composed of a sum of plane surface 
waves (henceforth denoted by SW3). Eq. H5.30j) tells us that the amplitudes 
dG^ of the SW3 in 73 decrease exponentially as the focal distance (i.e., /i+x^) 
increases, so that the SW3 make themselves felt all the less the farther the 
source is (in the vertical direction) from the ground. Note however, that the 
wavevectors associated with the SW3 are not identical to those associated 
with the SW2 because ki spans an interval in 73 that is different from the 
one in In addition, we note that to each horizontal wavenumber ki in 
the intervals ] — 00,— k^] and [k'^jOol, correspond a SSW in Oi and a SW3 
in f^Q . 

The main conclusion of this discussion is that for focal distances of the 
source that are sufficiently large, the diffracted field in the half space is 
essentially given by and is expressed by a sum of BW. This corresponds 
more or less to the situation in the quasi-lD analysis of the forward-scattering 
problem, but, as we shall see further on, it is, by no means, a valid picture 
of the response of the configuration when the focal distance of the source is 
not large. 

5.3 Amplitudes of the SBWl 

Henceforth, we restrict our attention to the field in the soft layer, and, in 
particular, to the three individual types of standing waves (SBWl, SBW2, 
SSW) of which it is composed. Here, we focus on a generic SBWl and note 
that its amplitude dG\ is the product of three factors: the factor S{u)) asso- 
ciated with the spectrum of the incident pulse, a geometric factor associated 
with the location of the source (whose influence was already discussed), and 
a so-called interference factor F^dki. We first discuss F^dki and then close 
the discussion with some remarks on S{u)). 

We rewrite ll as 



7|(x,a;) 



L 



F, 



{ki,uj) cos {ki{xi - xD) COS {K^iio)x2) e-^-^^HC^ 



IT 



(5.31) 



and are therefore interested in 



7^1 



\ {ki,uj)dki 



dk 



iK^{uj) cos(A'i(w)/i) + ^Kl{uj) sm[Kl{oj)h) 

ki e [0, 



(5.32) 



We make the change of variables 



ri = k'^h = h 




(5.33) 
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and examine Fl in the interval C G [0, 1]: 



lip cosycpr]) + v(psm.{(pr]) 



wherein: 



^^ = ^,7 = ^ = ^,^^ = , = a/^^ . (5.35) 

Note that 7 > 1 and f < 1 due to previous assumptions. Since rj and Q are 
real, the denominator in Fl cannot vanish; however it does attain minima 
for certain values of these parameters. 

Let us consider C, to be constant and inquire for what values of rj 

\Fl{Q,ri)\-'^ = ij^cos^icpT]) +v^(l)'^sm^{(l)r]) , (5.36) 

attains its minima. A necessary condition is: 

dv {\FHC,vT^) =0 = 4> {v'<t>' - V-') sm(2</.r?) . (5.37) 

There are three possibilities, the first one of which is = 0, but this implies 
C = 7 > 1 which is in contradiction with the fact ( must lie in [0, 1]. The 
second possibility is that ip = vcp] we will re-consider this case further on. 
The third possibility is sm{2(j)r]) = whence cpr] = n7r/2 ; n = 0, 1, .... To 
determine for what values of n these roots correspond to actual minima of 
1-^1 (Ci^) I ~^ we must have 

(|Fi(C,r?)|-2) l^^^^ = 2ct>^ (^;V' - cos(n7r) > . (5.38) 

This condition gives rise to two types of solutions depending on the sign of 
(f)^ — if)'^ . The first type, which we call even body wave solutions (designated 
by the superscript Be) is: 

V = Vm=—r ; m = 0,1,2,... when t;(/> > -0 • (5.39) 
(P 

The second type, which we call odd body wave solutions (designated by the 
superscript Bo) is: 

rj = r]E° = ^-^"^ ^''^ ; m = 0, 1, 2, ... when t;</>< i/; . (5.40) 

Let be the value of C, for which v(j) = ip. We find 



1 — f 



or 
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from which it follows that < 1, this meaning that the second possibility 
(i.e., V(j) = Tp) is not contradictory with the constraint C £ [0, 1]. 
Thus, the three types of solutions leading to minima of are: 



ForC>C'': ^ = ^m =^ ; m = 0,1,2,.... (5.43) 
For C = C^ ■■ all , (5.44) 

ForC<C'': V = V^"= ^ ' ; m = 0,l,2,.. . (5.45) 

The meaning of all this is that \fI\ has regularly-spaced (in terms of rj) 
maxima for all values of C, which is another way of saying that \Fi\ is a 
periodic function of rj for all (. The period of this function is ir/cp (even when 
C = C^: because a constant is a periodic function with arbitrary period). 
However, the function takes different forms in the three cases f H5.4Hjl - (|5.45p ). 
In fact, 

(i) for C > C^' has maxima equal to V ""^ = at = mir/cj) and 
minima equal to (u0)~"^ at rj = {2m + l)7r/2(/), 

(ii) for C = '1-^1 I ^ constant equal to tl^~^ = {v(j))~^ at all rj, 

(iii) for < '1-^1 I has minima equal to at rj = mir/cj) and maxima 
equal to {v(j))~^ at rj = (2m + l)7r/2(/). 

A numerical example will help to give a measure of the relative impor- 
tance of these three types of solutions. Recall that: 

^ (cW 

SO that 



(5.46) 



IS 



= V = • (^-^^^ 

Let us choose parameters that might be pertinent in the context of topics 
(a) and (b): c° =1000m/s, p° =1500kg/m3, =100m/s, p° =1000kg/m3, 
for which v = 0.67 x lO'^ and v-/ = 0.67 x 10~\ whence = 0.995. Thus, 
\Fi \ takes the form of the type (iii) function in most of the interval [0, 1], in 
fact in < C < 0.995. In particular, for body waves whose wavevectors are 
nearly- vertical (i.e., < C << 1), the maximum of 

1 ^ _ 

which, in the present numerical example, is equal to 15. 

The lowest frequency [v = u)/2tt) for which this value is attained (ob- 
tained from 1] = 7r/20 ~ is 

.= 1^, (5.49) 
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and is often called either the 'fundamental Haskell resonance frequency' 
[30] or the 'one- dimensional resonance frequency' [HSl, QB] of the soft soil 
layer/hard substratum configuration. However, a sinusoidal response func- 
tion of the type is not consistent with resonant response (which is infinite 
at the resonance frequencies in the absence of a dissipation mechanism) that 
would arise, for instance, in the context of excitation of some sort of struc- 
tural mode; in fact, this sinuoidal response results from interference of waves, 
which is the reason why we termed the 'interference factor'. Thus, it is 
improper to employ the term 'resonances' |54| . |53| in connection with 
body wave response (embodied in I^) of the configuration. 

To conclude this discussion, we now consider the spectral factor S{uj). It 
is obvious that if S{uj) = S{r]c^/h) is significantly large near the frequencies 
T] = (2m + l)7r/2(/) ; m = 0, 1, ... at which is large, then the product of 
these two functions, embodied in /J will be large at these frequencies. In 
particular, if S{uj) = S{r]c^ /h) is maximal near the low frequency rj = vr/2</>, 
then the response will be large over a large range of horizontal wavenumbers 
due to the contribution of the m = maximum of the interference factor 
F/. This has been noted repeatedly in the past (HOI, [^3]; and termed 
a 'resonant response' (as mentioned above), which it is not because t] = 
{2m + l)7r/2(/) ; m = 0,1, ... are not resonance frequencies. 

5.4 Amplitudes of the SBW2 

We again restrict our attention to the field in the soft layer, and, in 
particular, to the SBW2 component. We note that the amplitude dG^ of 
the generic SBW2 is the product of three factors: the factor S{u}) associated 
with the spectrum of the incident pulse, a geometric factor associated with 
the location of the source (whose influence was already discussed), and a 
so-called interference factor F2dki. Here we discuss the product Fj^dfci with 
S{uj) in order to evaluate the contribution of generic SBW2 to the overall 
response in the layer and on the ground. 

We rewrite I2 as 



h (x, t^) 
-S{lo) 




F2{ki,uj) cos {ki{xi — xD) cos [K2{uj)x2) e' 



TT 



(5.50) 



wherein 



F2 {ki,uj)dki 



dki 



— K. 



ki G [k^,i 



M] . (5.51) 
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We make the same change of variables as in the previous section, with the 

additional definition 

d := VC^- f , (5.52) 

and examine F2 in the interval ( G [1,7]: 

dc 



' cos((/)r/) + v(t) sm{(f)r]) 



(5.53) 



Contrary to the previous case, here the denominator in F2 can vanish for 
real rj and C, i.e., 

—9cos{(j)ri) + v(j)sm{(j)7]) = , (5.54) 

this being none other than the dispersion relation of Love modes. The roots 
of this relation are: 



1 



V 



arctan — - 


+ rrnr 


\V(I)J 





; m = 0,l,2 , 



(5.55) 



wherein the arctan function is defined in [— 7r/2,7r/2] and can be expressed 
either by the series 



arctan y = y + y~^(— 1 



.1 y 



21+1 



1=1 



21 + 1 



(5.56) 



or by the series 



(2^+1) 



arctan y = > — 1 — ; \ y > I 

^ 9 '21 + 1 



1=0 



(5.57) 



It is easily shown that 9 = vcj) when 



C = C'^ := A l + (72-l) 



(5.58) 



so that C'^ > 1, as it should be for the constraint C G [1,7] to be satisfied. 
Thus, three types of solutions lead to a zero in the denominator of 



For C < C : r] = 7], 



e__i 

v(h ?> \v 



+ 



T (4m + l)7r 
For C = C'' : n = vL = 



(5.59) 
(5.60) 



For C > C '■ V = "H: 



_ (2m + l)7r 1 



+ ■ 



v4> I f 
T ^ 3 



+ .... 



, (5.61) 
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and correspond to the existence of three types (even, neutral, odd) of Love 
modes whose eigenfrequencies are rj^, rj^, rj^ respectively. 

This means that | F2 | has regularly-spaced (in terms of rj) maxima for all 
values of C, which is another way of saying that \fI \ is a periodic function of 
r] for all C- The period of this function is ir/cj) (even when C = because a 
constant is a periodic function with arbitrary period). However, the function 
takes different forms in the three cases H5.59jl - (|5.61|l . In fact, 

(i) for C < C'^: 1-^2^1 has maxima equal to 00 at r/ = rj^, 

(ii) for C = -1-^2^1 has maxima equal to 00 at = 7/4; 

(iii) for C > :|^2^| has minima equal to 00 at = t]^. 

A numerical example will help give a measure of the relative importance 
of these three types of solutions. Let us again choose: =1000m/s, = 
1500kg/m^ =100m/s, =1000kg/m3, for which 7 = 10, -u = 0.67 x lO'^ 
and vj = 0.67x 10-\ whence = 1.0044. Thus, |F2^| takes the form of the 
type (iii) function for most of the interval [1,7], in fact in 1.0044 < 1^ < 10. 
A few remarks are in order. 

(1) contrary to what may be inferred from works such as [30], [131, i^ - 
|26| . the individual Love modes do not have the structure of surface waves 
in the layer (and, therefore, on the ground) since the SBW2 are actually 
standing body waves; the only feature they share with surface waves (i.e., 
the SW that coexist in the hard substratum when Love modes are excited) 
is their phase velocity 

= ^ , (5.62) 

wherein it can be noted that due to the fact that ( £ [1, 7], 

< c° , (5.63) 

which means that the phase velocity of Love modes (shared by the SBW2 in 
the layer and the SW in the hard substratum) is less than the phase velocity 
of body waves in the the hard substratum, 

(2) contrary to the what occurs in connection with the SBWl, the exci- 
tation of Love modes is indeed a resonant process, because Love modes are 
actually structural modes of the soft layer/hard substratum configuration 
and because the response associated with each of these modes is infinite at 
resonance in the absence of dissipation in both of the media of the configu- 
ration (this response can be large, but finite, when dissipation is present), 

(3) the resonant frequencies of the Love modes are not identical to the 
frequencies at which the SBWl attain their maxima; for instance, the differ- 
ence of these frequencies, for the m-th prevalent odd-type SBWl and SBW2, 
is: 

1 

3 



Bo _ Lo 
Im im 



V(j) 1 f V(j)^^ 



+ 



(5.64) 
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which means that the frequency of occurrence of the maxima of the m- 
th order SBWl is higher than (although it can be close to) that of the 
corresponding SBW2 (note that the difference in H5.64j) does not depend on 
m). 

To conclude this discussion, we again consider the spectral factor S{u)). It 
is obvious that if S{uj) = S{r]c^/h) is significantly large near the frequencies 
7/4° at which F2 is large (infinite if no dissipation is present), then the product 
of these two functions, embodied in I2, will be large at these frequencies. In 
particular, if S{i]c^ /h) is maximal near the low frequency t^q", then the 
response will be large over a large range of horizontal wavenumbers. If 
S{u}) = S{r]c^/h) is maximal near the low frequency 1]^°, and tjq" is not 
too far from tjq", then the global response can be even larger due to the 
cumulative contribution of both the SBWl and SBW2. 



5.5 Amplitudes of the SSW 

We are again concerned with the field in the soft layer, and, in particular, 
with its SSW component. We note that the amplitude of the generic 
SSW is the product of three factors: S{uj) which is associated with the spec- 
trum of the incident pulse, a geometric factor associated with the location 
of the source (whose influence was already discussed), and the interference 
factor F^dki. Here we discuss the product of F^dki with S{lo) in order to 
evaluate the contribution of generic SSW to the overall response in the layer 
and on the ground. 

We rewrite /| as 

4(x,u;) = 

r FUkuuj) cos {ki{xi-xl)) cosh (4{u)x2)e^>'^^''-'''^'^dki , 

(5.65) 

wherein 

dki 



Fl{ki,oj)dki 



K^{oj) COs]l{K\{L0)h) + ^^2(61-') sm\i{K\{ijo)h) 

ki G [A:\oo[ . (5.66) 



We make the same change of variables as in the previous two sections, with 
the additional definition 

X := , (5.67) 

and examine F^ for C in the interval [7,00]: 

FhC,v)dC = w . w s ■ (5.68) 

U cosh(x?7) + vx smh(x?7j 
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Since x ^ ^ ^oi C £ [7,00], and > 0, sinh(x?7) > and cosh(x??) > 
for C G [7,00], which means that the denominator in the previous formula 
cannot vanish for real r] and It can however exhibit minima for ^ G [7, 00]. 

Let us consider C to be constant and inquire for what values of ry the 
denominator (-^3^) = Ocosh^xv) + t'x sinh(x??) has minima. This requires 
that 

drj {Fi iCv)) =x[d sinh(x??) + vx cosh(xr?)] = . (5.69) 

But [ ] 7^ except for x = 0; i-^-, for ^ = 7 and Vr/. When x = we find 
(F^) ^ = 9, and from the fact that sinh(x??) > and ^cosh(x??) > for 
C G [7, 00], we conclude that (Fg^) ^ > 9. This means that C = 7 corresponds 
to the location of a minimum of (F^) ^ and this holds for all rj. 

Thus, F^ is a monotonically-decreasing function of C for all C £]7iC>o[ 
and attains its maximum equal to 9"^ = 1/ \/7^ — 1 at C = 7 for all t]. 

To get an idea of the magnitude of this function, notably in relation to F^, 
we again consider the numerical example: c" =1000m/s, =1500kg/m^, 
=100m/s, p° =1000kg/m3, for which 7 = 10, whence maxlFg^l < 0.1005 
which is much less than max | = 15 for the same set of parameters. 

Since the maximum of F^ is attained at all frequencies (i.e., for all rj), 
the spectrum function S{u)) does not influence the relative contribution of /g 
to the ground response. Thus, to conclude this discussion, we can say that 
the SSW contribute relatively little to the ground response in comparison to 
the SBWl and SBW2, except perhaps at frequencies close to the minima of 
the functions F^ and F2. 

6 Total frequency domain contributions of the SBWl, 
SBW2, SSW as embodied in the cumulative fre- 
quency response functions /|, I2 and for elastic 
and viscoelastic layers 

Although the theoretical analysis carried out in the sect. 5 may be useful 
for underlining the role played by the different types of body and surface 
waves that appear in the fields in the layer and substratum, it does not 
resolve the practical problem of the actual evaluation of the integrals /| , I2 
and I3. Another drawback of this analysis is that it is restricted to the case 
in which the layer is elastic, but the conclusions that were drawn for the 
elastic layer case should not be radically different for the case of a weakly- 
or moderately-viscoelastic layer. 

Consequently, we resorted to a purely numerical (i.e., Simpson integra- 
tion) approach for the evaluation of I2 and /| and of their sum to deter- 
mine the frequency- domain seismic response of the layer/substratum config- 
uration. Since physically-realistic configurations involve viscoelastic layers. 
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Figure 7.1: Moduli of source spectrum functions, i.e. |5(z^ = Lo'/27r)| versus 
i/(Hz), for I/O = 0.25Hz (left), 0.50Hz (middle), l.OHz (right). 

we evaluated these integrals and the total frequency response uixgjUj) under 
the assumption of viscoelastic layers. Once u{xg,uj) was computed, we de- 
termined the temporal signal u(xg,t), again by purely numerical means, via 

The weakness of the numerical approach is that it makes it difficult to 
discern the mechanisms underlying the observed response. To overcome this, 
we will give in sect. 8 a phenomenological analysis of the frequency-domain 
and time-domain responses which should facilitate the comprehension of the 
particular features of the temporal signals. 

7 Computational results 

7.1 Preliminaries 

In all except sect. 7.9 we take the density of the hard half space to 
be 2000kg/m'^. Contrary to what was assumed in the preceding theoretical 
analysis, we henceforth take into account the lossy nature of the soft layer. 
The quality factor is chosen equal to 30 in all the computations except 
in sect. 7.9 (recall that the hard half space is non-lossy, i.e., = oo). The 
seismic source is associated with the pseudo-Ricker impulse function given 
in H,S.4p whose spectrum is given in H,S.1()|) . Examples of these spectra (i.e., 
their moduli) are displayed in fig. 17.11 

Unless stated otherwise, the thickness of the layer h is taken to be 50m. 
This figure could just as well be 10km provided the wavelength and/or 
the wavespeed are adjusted so as to keep the ratio h/A^ = hu/c^ constant. 
This issue is discussed in more depth in the sect. 7.9. 

7.2 Comparison of the results of two methods for determin- 
ing the frequency domain response on the ground 

In order to be reasonably sure that the separation of variables technique 
employed herein gives valid results for a viscoelastic layer, we compared these 
results to those obtained by a finite element time domain viscoelastic code 
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Figure 7.2: Comparison of the frequency domain ground response (i.e., 
|«^(xg, V = uj/2'n)\ versus z^(Hz)) at x =(3000m,0m), for a shallow uq =0.5Hz 
source at x'^ =(0m,100m), in a Mexico City-like environment, i.e., = 
600m/s, = 1300kg/m3, = 60m/s. The full curve was obtained by the 
separation of variables technique described herein whereas the dashed curve 
was obtained by the finite element time domain technique described in 

inn. 



developed by one of the present authors (JPG) with C. Tsogka [HE], |3T] . 
The time domain responses obtained by this code were Fourier-transformed 
to get the corresponding frequency domain responses. An example of these 
results is given in fig. 17.21 



7.3 The cumulative contributions of the SBWl, SBW2 and 
SSW to the overall frequency domain ground response 

The discussion here centers on the transfer functions of frequency do- 
main ground response. In all that follows, the graph of the modulus of 
/J(xg, z^)/ii*(xg, z/) versus frequency u is designated by dots, the graph of 
the modulus of l\{'x.g,v)/u^{^g^i') versus v is designated by dashes, the 
graph of the modulus of (/2(xg,i/)^ + Ig (xg, z^))/u*(xg, zv) versus frequency 
V is designated by dot-dashes, and the graph of the modulus of the ground 
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Figure 7.3: Transfer functions of ground response in Nice-like environ- 
ment for various source locations and observation locations. Left: = 
(Om, 3000m), x = (100m, Om). Middle: x^ = (Om, 3000m), x = 
(3000m, Om). Right: x'^ = (Om, 100m), x = (3000m, Om). 

displacement u(xg , v) / {xg , u) versus frequency by a continuous line. 

To begin, consider a configuration thought to be representative of that 
in the central portion of the city of Nice (Prance) wherein = lOOOm/s, 

= 1800kg/m^, = 200m/s. We first place the source at a relatively- 
large depth of 3km on the X2 axis, i.e., x'^ = (Om, 3000m) and evaluate the 
moduli of the ground transfer functions relatively near the epicenter, i.e., 
X = (100m, Om) (left subfigure in fig. I7.3|l as well as relatively far from the 
epicenter, i.e., x = (3000m, Om) (middle subfigure in fig. 17. 3p . and then 
place the source at a relatively-small depth of 100m on the X2 axis, i.e., 
x^ = (Om, 100m) and evaluate the ground transfer functions relatively far 
from the epicenter, i.e., x = (3000m, Om) (right subfigure in fig. 17. 3p . 

It will be noticed that, in this and practically all subsequent results, the 
curve relative to [(I^x^,;^) -|- I-^Xg, j^))/tt*(xg, ;y)| is coincident with that 
relative to [/^(xg, J^)/ti*(xg, i^)| which means, as predicted by the analysis 
of the preceding section, that the contribution to overall ground response of 
the standing surface waves in the layer is negligible. Thus, we restrict the 
following discussion to the sole contribution of the standing bulk waves of the 
first (SBWl) and second kinds (SBW2). The left and middle panels in fig. 
l7.3l show that when the the focal depth is large the ground response is largely 
dominated by the contribution of the SBWl (i.e., by |/|(xg, i/)/n*(xg, z^)|), 
and, in fact, the SBW2 have no influence on the response beyond ~ IHz. 
However, the right panel in fig. 17.31 gives just the opposite result when 
the focal depth is small, since the ground response is largely dominated by 
the SBW2 (i.e., by |/2V5'(a;)|) and the SBWl have little influence beyond 
~ IHz. Another interesting feature of these results is that the total response 
curves have noticeably-different appearance when the source is deep or shal- 
low (notice that this appearance is qualitatively the same for small and large 
epicentral distances, assuming the same, large focal depths in the two cases). 

Next consider a somewhat softer environment than in Nice wherein = 
600m/s, pi = 1300kg/ m^, = 200m/s (fig. 17. 4p . We first place the source 
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Figure 7.4: Transfer functions of ground response in softer-than-Nice en- 
vironment for various source locations and at the fixed observation point 
X = (3000m, Om) . Left: = (Om, 3000m). Right: = (Om, 100m). 

(rather deep) at x* = (Om, 3000m) and evaluate the ground transfer func- 
tions relatively far from the epicenter, i.e., x = (3000m, Om) (left subfig- 
ure in fig. 17. 4|) . and then place the source at a relatively-small depth, i.e., 
X* = (Om, 100m), and again evaluate the ground transfer functions relatively 
far from the epicenter, i.e., x = (3000m, Om) (right subfigure in fig. 17. 4|) . We 
again observe that the ground response is dominated by the SBWl when the 
source is deep and by the SBW2 when the source is shallow. Also we notice 
that the appearance of the total response curve for a deep source is different 
from than of a shallow source. 

We next consider a Mexico-city like site (of course without the build- 
ings, contrary to the case in ^21) in which = 600m/s, = 1300kg/m^, 
= 60m/s (fig. I7.5|l . In all except the lower right hand panel we again ob- 
serve that the ground response is dominated by the SBWl when the source 
is deep and by the SBW2 when the source is shallow. The exceptional case 
is that of a shallow source and small epicentral distance, for which the con- 
tributions of the SBW2 and SBWl to the overall response are of comparable 
magnitude, especially near the first low frequency peak. A plausible cause 
of this behavior is the rather large contrast of body-wave velocities between 
the the layer and substratum, thus giving rise to a large contribution of 
the individual SBWl at the fundamental Haskell frequency (recall that this 
contribution is all the greater the greater the body wave velocity contrast). 

Next we consider a Mexico City-like environment with a somewhat harder 
substratum for which = 1500m/s, = 1300kg/m^, = 60m/s (fig. 17. 
For a deep source and large epicentral distance (left panel of the figure), the 
response is dominated, as usual, by the SBWl. When the source is shallow 
and the epicentral distance is large (right panel of the figure) we encounter 
a new kind of response characterized by contributions of the SBWl and 
SBW2 that are of comparable magnitude (this was obtained in the previous 
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Figure 7.5: Transfer functions of ground response in Mexico City-like en- 
vironment for various source locations and observation points. Upper-left: 
= (Om, 3000m), x = (3000m, Om). Upper-right: x'^ = (Om, 3000m), 
X = (100m, Om). Lower-left: x* = (Om, 100m), x = (3000m, Om). Lower- 
right: X* = (Om, 100m), x = (100m, Om). 
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Figure 7.6: Transfer functions of ground response in Mexico City-like with 
harder substratum environment for various source locations and at the fixed 
observation point x = (3000m, Om) . Left: = (Om, 3000m). Right: x* = 
(Om, 100m). 



figure for a shallow source and small epicentral distance. That this should 
occur even for a large epicentral distance is probably attributable to the 
fact that the body wave velocity contrast is very large (it was smaller in 
the configuration of the previous figure) which fact favorizes a substantial 
contribution of the SBWl (notably near the fundamental Haskell frequency), 
even when the distance between the source and observation point is large. 

The last result in this series concerns once again the Mexico City-like 
environment in which = 600m/s, = 1300kg/m^, = 60m/s (fig. 
17. 7p . We are now interested in evaluating the effect of changes in the layer 
thickness h for a shallow source and large epicentral distance. We observe 
in the figure that the response is dominated by the cumulative contribution 
of the SBW2 for all the layer thicknesses. Furthermore, the number and 
finesse of the resonance peaks in the interval [0, 2Hz] increases with h, the 
dominant peak always being the one associated with the resonant excitation 
of the first (lowest-frequency) Love mode and being located at a frequency 
that is all the lower the larger is h. 

7.4 Time records of various input pulses for different com- 
binations of source and observation point coordinates 

In the following, we exhibit (figs. I7.8I7.1()|) time records of the three 
pseudo-Ricker pulses having frequencies: vq = 0.25Hz, 0.5Hz, l.OHz (whose 
spectra were shown previously in fig. 17.11) . This is done for all combinations 
of the four coordinates (assuming xf = Om, X2 = Om): x| = 100m, j;| = 
3000m, xi = 100m, xi = 3000m. As expected, the pulses have the same 
shape for all source and observation point locations since the substratum 
was assumed to be an elastic (i.e., non-dispersive) medium; however their 
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Figure 7.7: Transfer functions of ground response in Mexico City-like envi- 
ronment for X* = (Om, 100m), x = (3000m, Om) and various layer thicknesses 
h. Upper-left: h = 20m. Upper- right: h = 40m. Lower-left: h = 60m. 
Lower- right: h = 90m. 
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Figure 7.8: Time records of the incident field (i.e., n*(x, t) versus t(s)) in 
the substratum (considered to fill all space and wherein (3^ = 600m/s) cor- 
responding to a 1^0 = 0.25Hz pulse for: x'^ = (Om, 3000m), x = (100m, Om) 
(upper left panel), x'^ = (Om, 3000m), x = (3000m, Om) (upper right panel), 
x^ = (Om,100m), x = (100m, Om) (lower left panel), and x'^ = (Om, 100m), 
X = (3000m, Om) (lower right panel). 
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Figure 7.9: Time records of the incident field in the substratum (considered 
to fill all space and wherein (3^ = 600m/s) corresponding to a fo = 0.5Hz 
pulse for: x'^ = (Om, 3000m), x = (100m, Om) (upper left panel), x** = 
(Om, 3000m), x = (3000m, Om) (upper right panel), x* = (Om, 100m), x = 
(100m, Om) (lower left panel), and x'^ = (Om, 100m), x = (3000m, Om) (lower 
right panel). 
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Figure 7.10: Time records of the incident field in the substratum (considered 
to fill all space and wherein = 600m/s) corresponding to a fo = l.OHz 
pulse for: x'^ = (Om, 3000m), x = (100m, Om) (upper left panel), x** = 
(Om, 3000m), x = (3000m, Om) (upper right panel), x" = (Om, 100m), x = 
(100m, Om) (lower left panel), and x'^ = (Om, 100m), x = (3000m, Om) (lower 
right panel). 
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Figure 7.11: Frequency and time domain representations of ground response 
in Nice-like environment for constant source-to-observation point distances; 
x'' = (Om, 3000m), x = (100m, Om) (solid line curves) and x* = (Om, 100m), 
x= (3000m, Om) (dotted line curves). Left: transfer functions. Right: time- 
domain responses (i.e., ii^(x, t) versus t{s)) to a fo = 0.5Hz pulse. 

maxima change as a function of these locations. Not unexpectedly, the 
largest pulses are those for which the source to observation point distances 
are the smallest. Of particular interest is the fact that the input pulse 
duration is approximately 2/fo, which corresponds to ~8s for the 0.25Hz 
pulse, ~4s for the 0.5Hz pulse, and ~2s for the l.OHz pulse. As will be 
seen hereafter, the response to these pulses in the layered configuration is 
generally of much longer duration. 

7.5 Comparison of frequency and time domain responses for 
constant source-to-observation point distances 

Again, we begin by a configuration thought to be representative of that 
in the central portion of the city of Nice (France) wherein c° = lOOOm/s, 

= 1800kg/m^, = 200m/s. Two constant source-to-observation point 
distance situations are considered: a) x* = (Om, 3000m), x = (100m, Om) 
(solid line curves in fig. 17. lip , and b) x^ = (Om, 100m), x = (3000m, Om) 
(dotted line curves in fig. 17.1111 . We notice in the left panel of fig. 17. Ill that 
the first bump of the transfer function occurs at a lower frequency when the 
source is near to the layer than when it is far from the layer, which fact 
suggests that the lower-frequency bump is due to the (resonant) excitation 
of the fundamental Love mode (SBW2) whereas the higher-frequency peak is 
associated with the first (non-resonant) interference (SBWl) maximum. The 
same remarks apply to the higher-order bumps. Moreover, the value of the 
transfer function at the first couple of bumps is much larger due to Love mode 
excitation than to constructive interference effects, and since the widths of 
these two bumps are approximately the same, the finesse (height/width ratio) 
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Figure 7.12: Time domain ground response in Nice-like environment for 
constant source-to-observation point distances and input pulse with vq near 
lowest-frequency maximum of transfer function. Left: = (Om, 3000m), 
X = (100m, Om), uq = l.OHz. Right: x^ = (Om, 100m), x = (3000m, Om), 
uq = 0.9Hz. 

of the Love mode peak is larger than that of the interference peak. The 
translation of this into the time domain is that the signal associated mainly 
with the fundamental Love mode resonance is more intense and of longer 
duration than the signal associated mainly with the fundamental interference 
bump. 

This last remark should be tempered by consideration of the spectrum of 
the input pulse, since the transfer functions do not take this spectrum into 
account whereas the temporal signals do. Thus, when the location of the 
maximum of the spectrum of the input pulse is closer to the location of the 
maximum of the transfer function, the time-domain response is larger, as seen 
in fig. 17.121 this being true for signals that are essentially due both to Love 
resonances or to constructive interference effects (note that the location of the 
pulse maxima were adjusted so as to be close to the locations of the transfer 
function fundamental peaks). Actually, this figure reveals the existence of a 
beating phenomenon in the ground response temporal signal for a source near 
the layer, which is probably due to the combined (amplitude modulation) 
effects of the fundamental Love mode peak and the fundamental interference 
peak. This issue will be discussed in more depth in the next section. 

Next consider the somewhat softer-than-in-Nice environment wherein 
c° = 600m/s, = 1300kg/m3, = 200m/s (fig. 171!^ . Again, the two 
source/observation point couples are: x^ = (Om, 3000m), x = (3000m, Om) 
(solid line curve in fig. 17.131 and x* = (Om, 100m), x = (3000m, Om) (dotted 
line curve in fig. I7.13|l . All that was said in the previous example concerning 
the transfer functions holds in the present case. Likewise, the repercussions 
on the temporal signals are the same as in the previous case (short duration 
signal for a remote source and relatively-long signal for a near source. 
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Figure 7.13: Temporal record of ground response in softer-than-Nice environ- 
ment for fixed source-to-observation point distance. The full curves in left 
(transfer function) and right (time domain response to fo = 0.5Hz pulse) 
panels refer to x'^ = (Om, 3000m), = (Om, 100m). The dotted curves in 
left (transfer function) and right (time domain response to = 0.5Hz pulse) 
panels refer to x* = (Om, 100m), x = (3000m, Om). 

We next consider a Mexico-city like site (again, without the buildings) in 
which = 600m/s, = 1300kg/m'^, = 60m/s (fig. I7.14p . Again, the two 
source/observation point couples are: x^ = (Om, 3000m), x = (3000m, Om) 
(solid line curve in fig. I7.14|) and x* = (Om, 100m), x = (3000m, Om) (dotted 
line curve in fig. I7.14() . All that was said in the two previous examples 
concerning the transfer functions holds in the present case. Likewise, the 
repercussions on the temporal signals are the same as in the previous two 
cases (short duration signal for a remote source and relatively-long signal for 
a near source. 

7.6 Time domain responses for very large and very small 
source-to-observation point distances 

We again consider a Mexico-city like site at which c° = 600m/s, = 
1300kg/m3, = 60m/s (fig. ITTTK)! . Although in fig. O upper left and 
lower right panels) we observed that the transfer functions for very large 
and very small source-to-observation point distances are qualitatively very 
similar, we are somewhat surprized to find that the two corresponding signals 
in fig. 17.151 are so qualitatively similar, due to the fact that the transfer 
function corresponding to the left panel in fig. I7.5l is dominated by the SBWl 
contribution, whereas the transfer function corresponding to the right hand 
panel in fig. 17.51 has strong contributions from both the SBWl and SBW2. 
The clue to this unexpected result resides in the spectrum of the i^o = 0.5Hz 
input pulse (see middle panel of fig. 17. 1|) . since the maximum of the latter 
is around vq = 0.6Hz and this frequency is both far-removed form the peaks 
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Figure 7.14: Transfer functions (left panel) and temporal records for vq = 
0.5Hz input pulse (right panel) of ground response in Mexico City-like 
environment for various source locations and observation points: x'^ = 
(Om, 3000m), x = (100m, Om) (solid line curves), x^ = (Om,100m), x = 
(3000m, Om) (dotted line curves). 
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Figure 7.15: Temporal records of ground response in Mexico City-like en- 
vironment for various very large (left panel: x* = (Om, 3000m), x = 
(3000m, Om)) and very small (right panel: x* = (Om, 100m), x = 
(100m, Om)) source-to-observation point distances). 
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Figure 7.16: Temporal records of ground response for = 0.3Hz (left panel), 
vq = 0.5Hz (middle panel), fo = 0.9Hz (right panel) input pulses in Mexico 
City-like environment for = (Om, 100m), x = (3000m, Om). 

of the transfer functions and characterized by a predominant contribution of 
the SBWl (the latter fact providing an explanantion of the relatively-short 
duration of the response signals in fig. 17.151 (note that the intensity of the 
very close source-to-observation point signal is much larger than that of the 
other signal, as it should be). 

7.7 Time domain responses for different input pulses and 
fixed source and observation point coordinates 

It is of considerable interest to ascertain to what extent the spectrum of 
the source affects the ground response ^Hl, notably when the source is near 
the layer. In fig. I7.1()l (which again applies to the Mexico City like environ- 
ment), we see, as expected, that when the spectrum of the input pulse is such 
as to overlap substantially the frequency band covered by the fundamental 
Love mode peak (left and middle panels in fig. I7.16|l . the duration of the 
signal is large. When the spectrum of the input pulse is such as to over- 
lap substantially the frequency band covered by the second Love mode peak 
(right panel in fig. I7.16|l . the duration of the signal is smaller than in the 
previous case due to the smaller finesse of the second Love mode resonance 
peak. In all three cases, we observe some beating, presumably due to the 
proximity of the interference peak and Love resonance peak in the transfer 
function. 

7.8 Time domain responses for various layer thicknesses and 
input pulses 

A more thorough study of the infiuence of the input spectrum must take 
into account variations of the layer thickness. This is done (again in 
the Mexico City-like environment and for a source that is 100m below the 
ground) in fig. 17. 171 for a 20m thick layer, in fig. 17. 181 for a 40m thick layer, 
in fig. 17. 191 for a 60m thick layer, and in fig. 17.201 for a 90m thick layer. We 
observe quite different responses, varying from a short pulse quite similar 
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Figure 7.17: Temporal records of ground response for uq = 0.25Hz (left 
panel), z/q = 0.5Hz (middle panel), vq = l.OHz (right panel) input pulses 
in Mexico City-like environment with layer thickness h = 20m for = 
(Om, 100m), X = (3000m, Om). 




Figure 7.18: Temporal records of ground response for = O.lHz (left panel), 
z/Q = 0.25Hz (middle panel), = 0.5Hz (right panel) input pulses in Mexico 
City-like environment with layer thickness h = 40m for x^ = (Om, 100m), 
x= (3000m, Om). 




Figure 7.19: Temporal records of ground response for vq = O.lHz (left panel), 
z/Q = 0.25Hz (middle panel), = 0.5Hz (right panel) input pulses in Mexico 
City-like environment with layer thickness h = 60m for x* = (Om, 100m), 
x= (3000m, Om). 
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Figure 7.20: Temporal records of ground response for vq = O.lHz (left panel), 
vq = 0.25Hz (middle panel), vq = 0.5Hz (right panel) input pulses in Mexico 
City-like environment with layer thickness h = 90m for = (Om, 100m), 
X = (3000m, Om). 

to the input pulse (for the thinnest layer and the lowest frequency input 
pulse) to a very long duration pulse (as much as 200s as compared to the 
the 4s duration of the input pulse) with pronounced beating (for the thickest 
layer and a medium frequency input pulse). Note that the 90m layer also 
corresponds to the case in which the source is closest to the layer (10m from 
the bottom face of the layer), which may also be a factor contributing to the 
pronounced anomalous character of the response in this configuration. 

Finally, we consider a Mexico City-like environment with a somewhat 
harder substratum for which = 1500m/s, = 1300kg/m^, = 60m/s 
(fig. 17.2111 . Two constant source-to-observation point distance situations are 
again considered: a) x'^ = (Om, 3000m), x = (100m, Om) (solid line curves 
in fig. [7^ . and b) x'' = (Om, 100m), x = (3000m, Om) (dotted line curves 
in fig. 17.2111 . For a deep source and large epicentral distance, the frequency 
response (in the left hand panel of the figure) is dominated, as usual, by the 
fundamental interference peak, and this results in a relatively-long duration 
time domain signal (solid line curve in the right hand panel of the figure) 
in response to a z^o = 0.5Hz input pulse. When the source is shallow and 
the epicentral distance is large the frequency response is dominated, unsur- 
prizingly, by the fundamental Love mode resonance peak, and this gives rise 
to a somewhat longer-duration signal (dotted line curve in the right hand 
panel of the figure) in response to the uq = 0.5Hz input pulse. This example 
indicates that it may be difficult to distinguish between the contributions of 
the SBWl and SBW2 when the contrast of the material properties between 
the layer and the substratum is very large. 

7.9 Regional path effects 

We now examine the manner in which a wave radiated from a source 
located underneath, but close to, the lower crustal boundary propagates 
over long distances. 

An example (real) of such motion, relative to the 11-3-02 (shallow) seis- 
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Figure 7.21: Transfer functions (left panel) and temporal records (right 
panel) of ground response to a fo = 0.5Hz input pulse in Mexico City- 
like with harder substratum environment for various source locations and 
observation points = (Om, 3000m), x = (100m, Om) (solid line curves), 
x^ = (Om,100m), x = (3000m, Om). 

mic Denali (Alaska) event recorded at a free-field ground location (i.e., 1.5km 
from the buildings of the city of Anchorage) at an epicentral distance of more 
than 275km is given in fig. 10 of [H]. A remarkable feature of this motion is 
its long duration of over 125s. We will show that this is possible with our 
simple model. 

We constructed a hopefully-plausible crustal model starting with the 
parameters of a thin layer, softer-than-Nice like configuration (for which 
po = 2000kg/m3, pi = 1300kg/m3, c° = 600m/s, = 200m/s, Q° = oo, 

= 30, /i = 80m), and by assuming conservation of such quantities as 
k^h, p^/p^, etc. in going to a much thicker layer. Let us suppose that 
we have two configurations, one of which is thin-layered and known (con- 
figuration with subscript 1), and the other is thick-layered and unknown 
(configuration with subscript 2). The layer in the known configuration is 
relatively soft and lossy, whereas it is relatively hard (although always softer 
than the substratum) in the unknown configuration. Since harder media 
are usually less lossy, we assume rather arbitrarily that the Q of the layer 
in the unknown configuration is 20 times larger than the Q in the known 
configuration, while the Q's of the substratum remain infinite in both con- 
figurations. Thus, we have: Q2 = Qi = oo and = 20Q\ = 600. Conser- 
vation of k^h means /cj/ii = A;2/i2, whence C2 = c}f2/i2/j^i^i, or, if we choose 
vi = IHz, f2 = 0.08Hz, and /i2 = lO/cm, then C2 = 2000m/s. Conserva- 
tion of p^/p^ means Pi/p\ = P2/P21 ^o that if we choose P2 = 2600kg/m'^ 
(close to the density of granite), then P2 = 1690kg/m^. We would also 
like to conserve wavespeed proportions, i.e., Ci/c\ = 02/0^, but this turns 
out to give wavespeeds in that are much larger than those for gran- 
ite (3200m/s) for the choice C2 = 2000m/s, so that we arbitrarily chose 
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Figure 7.22: Ground motion (displacement) at large epicentral distance 
(xi — xf = 300km) in response to a pseudo-Ricker pulse line source un- 
derneath, and close to (focal depth x| = 12km), the lower boundary of 
a thick {h = 10km), fairly hard, crust overlying a granite-like substratum 
(po = 2600kg/m3, pi = IGOOkg/m^, c° = 3000m/s, = 2000m/s, Q° = oo, 
= 600, h = 80m). Transfer functions, with same notations as in fig. 17.31 
(upper left panel) and time histories for: & uq = 0.05Hz input pulse (upper 
right panel), a f o = O.lHz input pulse (lower left panel), a f o = 0.2Hz input 
pulse (lower right panel). 

C2 = 3000m/s (i.e., close to the wavespeed in granite). Finally we chose 
to conserve the relative distance of the source to the lower boundary of the 
layer, i.e., (x| — h)/x2, which gives rE| = 12km. The results for this new 
(crustal) configuration excited by the usual pseudo-Ricker pulse line sources 
are given in fig. 17.221 wherein it can be seen that the ground response far 
from the epicenter (300km): 1) is dominated by the excitation of Love modes 
(notably the fundamental), 2) takes the form of a pulse which has a shape 
not very different from that in fig. 10 of (H] and is of approximately 100s du- 
ration, governed essentially by the fundamental Love resonance peak. Thus, 
our theoretical model shows that it is quite possible for a source underneath, 
and relatively close to, a fairly thick, fairly hard crust overlying a very hard 
substratum to give rise to a rather long-duration pulse even at large epicen- 
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tral distances. What becomes of this pulse when a city is located at this 
large lateral distance from the source constitutes an important, and as yet 
not fully-elucidated, question (this meaning, that although studies such as 
[HI are designed to take into account all that occurs between the distant 
source and the observation point in the basin, the results that are offered are 
entirely of numerical nature and therefore do not provide an explanation of 
the underlying physical processes). 

8 Phenomenological model of the time-domain 
ground response 

In sect. 6 we mentioned the difficulties of obtaining closed- form expres- 
sions of the integrals ll, and I3 and therefore of the Fourier integral 
H3.35jl accounting for the time-domain ground response. Nevertheless, the 
many numerical examples in sect. 7 of the frequency- domain ground re- 
sponse all seem to have common features which we shall attempt to describe 
in this section in phenomenological manner. Moreover, this approach will 
be shown to provide a simple means of understanding the origin of the main 
features of the time domain response. 

The principal features of the ground transfer function |u(xg, u;)/ii*(xg, u;)| 
appeared to be due to interference and resonance causes and manifested 
themselves by a series of well-defined, regularly-spaced bumps. That the 
nature of these bumps be due either to interference or to resonance causes is 
not of primal importance at the present (phenomenological) level of analy- 
sis; the only aspects that interest us now are the relative widths and heights 
of the bumps (recall that, in general, the bumps associated with interfer- 
ences (SBWl) are broader and less intense than the corresponding bumps 
associated with Love mode resonances (SBW2)). 

We represent each of these bumps by a gaussian function of frequency, 
which we multiply by the spectrum S{oj) of the incident pulse and by an 
amplitude function ^(x^, x*, w) of x^, x^ and u; to take into account the fact 
(observed in the numerical results) that the different bumps of frequency- 
domain response indeed depend on these quantitites. Let Gi{uj) be the l-ih 
gaussian function of the form 

1 ("-"i)^ 

Gi{uj) = ^=e . (8.1) 

The bump connected with this function attains its maximum aA. uj = uji and 
its finesse is all the larger, the smaller is In fact (p. 319), Gi{uj) tends 
towards the Dirac delta distribution 6{lo — uJi) as £1 — > 0. 
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Thus, we represent the frequency-domain ground response by 

L 

« J]5HA(x„x^w)G^(a;) , (8.2) 
1=1 

wherein w^+i > uji and L may be a large positive integer. However, the latter 
will be taken equal to 4 due to the fact that we assume that the spectrum 
S{u;) of the input pulse is not too broad and centered at low frequencies (the 
case of interest in the applications considered herein). Introducing (|8.2|1 into 
l|3.35|l gives 

4 poo 

u^(xg,t) « 23f? V / S{uj)Ai{xg,x',Lo)Gi{u;)e-''^^dLJ . (8.3) 
1=1 -^0 

Although we don't know, nor assume, much about S and Ai, it seems rea- 
sonable to suppose that these functions are slowly-varying in comparison to 
Gi{uj) with respect to u; in the neighborhoods uji in which the gaussians are 
maximal. Consequently, we can make the approximation 

u\^g,t)^2^y^S{u;i)Aiixg,x',u;i) Gi{u:)e-'^'dio . (8.4) 
1=1 

By proceeding as in [HE] (p. 313) and making use of the identity [52] 
exp(— |^)(i^ = \/vr/2, we find 

/>oo 

/ Gi{uj)e-''^'duj = e-T*'-^-i* , (8.5) 
Jo 

so that 

4 

u\-Kg,t) « 25^^5(a;0^^(x„x^.;0e-^*'~*"'* . (8.6) 
(=1 

By representing S and Ai in polar form 

S(a;0 = |5(^z)|e*'^'(-') , A(x„ x^ c^,) = |^;(x„ x^ , 

(8.7) 

we get 

4 

u\^g,t) ^ 5? J]SKxg,x^u;z)e-T*'-N*-A(---^-')] = 
1=1 

4 

^B/(xg,x^a;i)e"~* cos (cj/t - /?/(xg, xs, cj;)) , (8.8) 
(=1 
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wherein 



e(x3,,x^u;0 := 2|5(a;^)M^(x3,x^u;^)|, 

/3^(xg,x^u;/) := ai{xg,x'' ,uji) + a{uJi) . (8.9) 

With this in hand, we now try to account for the main features of the nu- 
merical results pertaining to the time records of ground response. 

Assume that one of the bumps in the ground transfer function dominates 
all others. In the present paradigm, this signifies that one of the terms, say 
the m-th in the sum in H8.8|) dominates all the others. As observed in the 
numerical results, this term should account either for the fundamental Love 
mode resonance (m=l) or for the fundamental interference peak (m=2), the 
necessary condition for the Love peak to be dominant being that the source 
is close to the bottom boundary of the layer, and the necessary condition for 
the interference peak to be dominant being that the source is far from the 
bottom boundary of the layer. In either case, we have: 

U^{Xg,t) ^ Bm{Xg,X^,UJm)e~^''\os{uJmt - Pm{Xg,X.'',UJm)) , (8.10) 

which is indicative of the existence of a monochromatic wave (angular fre- 
quency LOm) whose amplitude Bm and phase /3m vary with the positions of the 
source and the observation point on the ground, and which is exponentially- 
attenuated with time. This attenuation is more or less pronounced, so that 
the duration of the signal is large if Em is small (i.e., the finesse of the cor- 
responding transfer function bump is large) as would occur for a Love mode 
resonance, and is relatively small if Em is relatively large (i.e., the finesse of 
the corresponding transfer function bump is relatively small) as would gen- 
erally occur for an interference peak. The same phenomenon is produced, 
although with less intensity due to the lowering of \Bm\, if the spectrum of 
the incident pulse is such as to favorize either the second Love mode res- 
onance or the second interference peak rather than the fundamental Love 
mode resonance or the fundamental interference peak. 

Consider a different type of situation in which two bumps in the ground 
transfer function dominate all others. In the present paradigm this means 
that two terms in (|8.8p . say the m-th and n-th, dominate all the others, i.e., 

U^(Xg,t) ^m(Xg,X^ u;m)e~^*^ cos {uJmt - ^m(Xg,X^ Wm)) + 

j-2 

Bn{xg,x\uJn)e~~ COS (a;„t - /3„(xg,x'',u;„)) . (8.11) 

The numerical results in the previous section show that this occurs only 
when the m-th term corresponds to a Love mode resonance and the n-th 
term to the interference peak nearest the Love mode resonance peak, so that 
n = m-|-l in the present numbering system. Moreover, the most pronounced 
phenomena were shown numerically to be produced when the fundamental 
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Love mode resonance and fundamental interference peak are involved, so 
that m = 1, n = 2 is the most interesting case. 
We can write (e.g., for m = 1, n = 2) 



n^(xg,t) 



Bie 4 
so that 



+ ^26- 



C1+C2 



+ [Bie 



626- 



C1-C2 



i.l2) 



uHxg,t)^ (Bie-"^'' +626-"} 



cos 



(Wi + UJ2)t (3l + 132 



cos 



(Wi - L02)t Pi - (32 



Bie 



£2 -,2 



sm 



{L0l+UJ2)t Pi +(32 



sm 



{UJI - UJ2)t Pi - P2 



§.13) 



An interesting case is: Bi ^ B2, £1 ^ £2, whence 



2B^e- 



cos 



{LJl+LJ2)t P1+P2 



COS 



(wi - u;2)t Pi - P2 



i.U) 



which is indicative of monochromatic sinusoidal motion of angular frequency 
'^^"^"^ , amplitude-modulated by an attenuated sinusoid of frequency | '^^"^2 1. 
In this case, the signal associated with the ground motion exhibits the beat- 
ing and attenuation observed in some of the computed results, with the 
duration depending, as in the previous case, on the finesse of the resonance 
and interference bumps. 

Actually, it is not necessary for Bi ^ B2, £1 ^ £2 in order to have beating 
in the signal, since although the motion associated with (j8.1,Sp is more com- 
plicated than that of (|8.14p . in both form of attenuated signal with 
more or less regular beating is present. The signal with irregular beating 
predicted by (|8.14p is the case most commonly observed in the computed 
results of the previous section. 

Thus, the phenomenological model accounts for all of the features of 
the computed time records: quasi-monochromatic attenuated signals (whose 
duration is governed by the finesse of the frequency domain bump) with- 
out beating when either the Love mode resonance or interference peak is 
involved, an attenuated signal with regular beating when the frequency do- 
main bumps are of comparable magnitude and finesse (the latter governing 
the duration of the signal), an attenuated signal with irregular beating when 
the magnitude and finesse of the fundamental Love and interference bumps 
are rather different. 
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9 Discussion 



We shall now attempt to provide answers to the questions raised in sect. 

1. 

The first question was: is it possible to obtain anomalous response with- 
out any lateral heterogeneity in the underground medium? The configuration 
studied herein was laterally- homogeneous. We have shown that ID response 
only accounts for interference effects (as embodied by the SBWl), but not 
for coupling to Love modes (as embodied by the SBW2) in the layer, which 
is particularly strong when the source is in the neighborhood of the lower 
boundary of the layer. Insofar as anomalous effects are essentially charac- 
terized by long duration and beating phenomena in the signals (e.g., curves 
in the middle and right panels of fig. I7.19|) . the answer to this question is 
negative as concerns ID response. However, when the contrast of material 
properties between the layer and substratum is very large, it is possible to 
obtain fairly-long duration signals (albeit without beating), which are es- 
sentially associated with ID response, even when the source is far from the 
lower boundary of the layer (solid curve in right panel of fig. 17.21(1 . More 
generally, i.e., when coupling to Love modes is achieved, the answer to the 
question is positive. 

The second question was: what is the relation of ID to 2D response and 
how adequate is it to model the general response of the configuration by its 
response to a (nearly) vertically-incident plane wave? We have shown that 
not only does the ID model not give rise to resonance phenomena, but that 
truly-resonant phenomena associated with the excitation of Love modes can 
only be described by a fully 2D (or 3D) model. For a source far from the lower 
boundary of the layer, the response is essentially due to the contributions of 
the SBWl (more or less equivalent to the ID response), but when the source 
is near this boundary, the waves (SBW2) not included in the ID model play 
a major role in the overall response in that they either overwhelm the ID 
response (long duration response without beating) or combine with the ID 
response to produce signals with long duration and beating. This finding 
should be taken into account in relation to studies (e.g., pi]) that attempt 
to predict seismic response of urban sites from ID type of analysis. 

The third question was: how does the focal distance of the source affect 
the response? The answer to this question was already provided in the previ- 
ous two paragraphs. However, it is opportune to reconsider this question in 
the light of topic (b) concerning the effects of underlying soil heterogeneities, 
lateral variations of the underlying soil layer, and built environment on seis- 
mic response in urban sites. One can show f|57p that a wave incident on a 
heterogeneous medium gives rise to a diffracted wave which can be considered 
to be radiated by induced sources (as opposed to the active source associated 
with the primary seismic disturbance) located within the medium. These 
induced sources can also appear on the boundary of the medium (especially 
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at endpoints, corners and irregularities of the boundary), so that the edges of 
a soft basin or the stress-free ground which includes the buildings overlying 
a homogeneous soft layer in a city-like site, can also constitute the locations 
of intense induced sources in response to an incident seismic wave. The fields 
radiated by all these induced sources can be represented in a manner identi- 
cal (provided the basic geometry of the configuration is similar) to that of the 
present work, so that much of what was written and found above, notably 
concerning the response to active sources located outside, and in the vicinity 
of the soft layer (and, by extension, to induced sources located within or on 
the boundaries of the soft layer, should apply to city-like sites built on soft 
layers or basins. The most important point (mentioned in references such as 
|12j . implicit in jSH], (2^, and proven herein as concerns active sources) is the 
following: the presence of these active or induced sources, located near or 
within the soft layer overlying a relatively-hard substratum, enables coupling 
to Love- type modes which are responsible for a part of the anomalous ground 
response observed in cities such as Mexico, notably motion characterized by 
long durations and beatings. 

The fourth question was: how does the epicentral distance affect the 
response? We have shown that the epicentral distance {\xi — xfl, for = 
(xi,0) is not as sensitive as other parameters (especially the focal distance 
|0 — ?/f I) as concerns its influence on duration and on the presence or absence 
of beatings in the ground motion (it should be mentioned that in the 
duration appears to be a linearly-increasing function of epicentral distance, 
but the slope of this function decreases when the crustal layer is softer). 
However, the epicentral distance (more generally: the distance of the source 
to the observation point) is a critical factor in determining the intensity of 
the signal on the ground (see, e.g., figs. 17. 14117. 15p . 

The fifth question was: how does the contrast of mechanical properties 
between the layer and the half space affect the ground response? There does 
not appear to exist a clear-cut answer to this question (see, however, [HHl 
in which it appears rather systematically that softer layers lead to longer 
durations and lower peak response for a given epicentral distance), since the 
dependence on the mechanical parameters is very much intermingled with 
that on the values of the geometrical and source parameters. On the whole, 
most of the answers to the previous four questions hold in a qualitative sense 
whatever the contrast of mechanical properties (see, e.g., the results herein 
for Nice, softer Nice, Mexico, Mexico with harder substratum), although 
there are some quantitative differences (e.g., the contrast influences consid- 
erably the intensity of the SBWl contribution). Naturally, the most remark- 
able features of the ground response of our simple configuration, which are 
due to the excitation of the fundamental Love mode, can only be observed if 
the layer is softer than the substratum in the sense of (|5.,Sjl and (|5.4|1 . These 
conditions are so broad and widespread (for city-like sites) as to render the 
'anomalous' response described herein a quite universal phenomenon. 
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The sixth question was: how does the thickness of the layer affect the 
response? In fig. I7.7l it was found that increasing the layer thickness increases 
the number of peaks, as well as the finesse of each of the latter, in a given 
range of frequencies of the transfer function. This has the effect of lowering 
the frequency of occurrence of the first peak so that the time-domain response 
will be largely conditioned by the spectrum of the input pulse, assuming the 
latter to be centered at a relatively high frequency. Thus, a low frequency 
pulse can produce substantially the same type of response for a thick layer 
as a relatively high frequency pulse in a thin layer. This point is important 
in connection with the topic of regional path effects mentioned in sect. 1. 

The seventh question was: how do the spectral characteristics of the in- 
cident pulse affect the response? The answer to this question can be found 
by comparing the three subfigures in any one of figs. 17.17117.201 Obviously, 
the spectrum of the incident pulse is a key factor (see sect. 8), which: a) if 
it overlaps either a constructive interference peak or Love mode peak, gives 
rise to attenuated, quasi monochromatic response, often of long duration 
(see figs. 17.17117.201 in which an example is given of a pulse having a dura- 
tion of 4s that gives rise to substantial ground response of 200s duration), 
b) if it overlaps both a constructive interference peak and Love mode peak, 
gives rise to attenuated, quasi monochromatic response with more or less 
regular beatings, c) if it doesn't overlap significantly either a constructive 
interference of Love mode peak, gives rise to a time domain response that 
can be qualitatively quite similar to the input signal (see, e.g., left panel of 
fig. I7.18j) . When the sources are induced, their spectra will be modified with 
respect to that of the spectrum of the primary active source due to diffrac- 
tion and dispersion, so that an a priori unfavorable situation for anomalous 
response from the point of view of the primary active source may turn out 
to be favorable from the point of view of the induced sources. 

10 Conclusion and perspectives 

This work originated in the observation that no satisfactory physical 
explanation has been given until now of anomalous seismic response in urban 
environments with soft layers or basins overlying a hard substratum. The 
principal reason for this knowledge gap probably lies in the complexity of the 
sites examined in previous (essentially-numerical) studies: 1) a homogeneous 
or multilayered basin of complicated form not including buildings (e.g., fH^ . 
[HI, EH]); 2) a homogeneous layer overlain by a periodic or non-periodic set of 
blocks or buildings (e.g., p2], (H^). The choice was therefore made herein to 
simplify as much as possible the characteristics of the site and sollicitation, 
while retaining as many as possible of their essential features. Thus, it was 
thought that: i) the problem had to be treated at least as a 2D one, ii) the 
sollicitation should not be a plane wave (for which coupling to Love modes 
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is impossible in the chosen configuration) but rather the wave radiated by 
a source which could be as simple as a line source ( this source eventually 
being able to mimick induced sources in more complicated configurations), 
iii) the soft component of the site could be a layer (rather than a basin) 
with flat, horizontal boundaries (i.e., flat rather than irregular ground, as 
rendered by the presence of buildings, flat interface with the substratum, 
rather than curved or irregular as for a basin or irregular layer). 

In spite of the simplicity of the model, obtaining an explanation of the 
principal features of the seismic response turned out not to be straightfor- 
ward. The temporal response took the form of a Fourier transform (with 
respect to frequency) of a frequency domain response function which itself is 
an integral with respect to the horizontal wavenumber component of plane 
body and heterogeneous waves. It was shown that the wavenumber integral 
splits quite naturally into three parts comprising either standing body waves 
in the layer (SBWl) and propagating body waves (BW) in the substratum, 
standing body waves in the layer (SBW2) and surface waves (SW2) in the 
substratum, standing surface waves in the layer (SSW) and surface waves 
in the substratum (SW3). It turns out that the amplitudes of the SW and 
SBW2 diminish exponentially as the vertical distance of the source to the 
lower boundary of the layer increases so as to make the contribution of the 
SBWl preponderant for sources with large focal depths. This fact provides 
an explanation of the relative success of the ID model (a remote source ra- 
diates a wave that has practically all the attributes of a plane wave when 
arriving on the layer) for remote sources, but also of the reason why the ID 
model is inappropriate for active (and, by extension, induced) sources that 
are close to (and, by extension, within) the layer, since the SBW2 waves 
are not included in this model. It was shown that each SBWl/BW pair is 
the principal ingredient of the ID model and that the maxima of the am- 
plitudes of the SBWl do not correspond to resonances, but are rather the 
result of the constructive interference of standing waves in the layer. Each 
SBW2/SW2 pair turned out to be a Love mode when the frequency satis- 
fies the Love mode dispersion relation. The amplitudes of each SSW/SW3 
pair were found to be negligible compared to those of the SBWl/BW and 
SBW2/SW2 pairs. 

The theoretical analysis reached its limits when the horizontal wavenum- 
ber integration was attempted. Thus, the integrals appearing in the frequency- 
domain response were carried out numerically and a parametric study was 
made of the cumulative contributions of the SBWl and SBW2. It was shown, 
as expected, that the SBWl give the preponderant contribution for remote 
sources, while both the SBWl and SBW2 cumulative contributions can be 
significant for nearby sources. The interference nature of the amplitudes of 
the individual SBWl was shown to be maintained in the frequency-domain 
cumulative response of these waves. The resonant nature of the amplitudes 
of the individual SBW2 was shown to be maintained in the frequency-domain 
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cumulative response of these waves. However, it was not possible to obtain 
mathematical expressions for the integrals of the frequency-domain responses 
appearing in the global time-domain ground response. 

We also showed that it possible for a source, underneath, and relatively 
close to, a fairly thick (10km), fairly hard crust overlying a very hard substra- 
tum, to give rise to a rather long-duration pulse even at large (e.g., 300km) 
epicentral distances, and that this finding is in agreement with what has 
been observed in connection with real earthquakes (see, e.g., (El, [HI)- We did 
not carry out an extensive analysis of this finding, nor address the issue of 
what becomes of this pulse when it enters an urban center located at large 
lateral distances from the source (as was done numerically in works such as 

ini,[ini,[ini,[ni). 

In the last section of this investigation, a phenomenological model was 
introduced based on the observation that the frequency- domain cumulative 
response components of both the SBWl and SBW2 appear as a series of 
regularly-spaced bumps which were modeled as gaussians. This enabled a 
closed form expression of the integral of the frequency- domain responses 
to be obtained which revealed and accounted for the type of time-domain 
response obtained by purely numerical means, notably, its attenuated, quasi 
monochromatic character, with regular or irregular beatings, governed by 
the finesse and relative position of the frequency domain bumps. This type 
of response has often been observed in earthquake-prone cities built on soft 
soil, so that it may be that some of the causal agents inherent in our simple 
model are operative in more complicated sites. 

A question that naturally arises is whether the type of analysis carried 
out herein can be extended to more realistic configurations in which induced 
sources are likely to play a major role. Our feeling is that this can be done 
provided some clever approximations are made in the expressions for the 
response of these configurations. 

Another question (alluded-to in one of the previous paragraphs) is that 
of regional path effects on global response in cities such as Mexico subject 
to earthquakes arising from laterally-remote sources. This very important 
theoretical issue, already considered in works such as ^^1, will have to be 
treated in more depth, first in the manner of the present contribution, to 
examine how the wave radiated by the source reaches the city site, what 
the nature of the waves are when they arrive in the city, and how they 
are converted therein into the form they have been observed to take (quasi 
Love or Rayleigh waves giving rise to high intensity, extremely long (even 
longer than what was found herein) duration ground motion, accompanied 
by beatings). 

Most of the extensions of the present work will have to be carried out 
first in the 2D, shear horizontal wave context in order to discern the essential 
issues. The extensions to the 2D- P/SV (as in e.g., jH]) and 3D (as in e.g., 
|59| ) cases with more general types of sources [H], ^H] are, of course, the 
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requisites for a full understanding of what happens when a seismic wave hits 
a realistic urban site. 
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